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ABSTRACT 


'The  chief  purpose  of  this  resort  is  to  provide  an  introduction 
to  pseudo-noise  modulation  functions  and  to  describe  those  proper¬ 
ties  that  make  their  use  in  radar  systems  desirable.  B^tause  of 
the  interests  of  these  laboratories,  characteristics  applicable  to 
fuzes  are  emphasized. 

The  applicability  and  usefulness  of  various  pseudo-noise  modu¬ 
lation  functions  in  radar  systems  is  discussed.  The  use  of  the  m- 
sequence  in  a  complete  radar  system  is  described,  and  an  extensive 
bibliography  of  other  applications  to  radar  systems  is  also  given. 

A  short  t-oction  containing  classified  information  on  jamming  of 
pseudo-noise  modulated  radars,  with  bibliography,  is  issued  as  a 
supplement  to  this  report. 

This  report  is  essentially  a  survey,  and  the  body  of  the  report 
therefore  contains  little  new  material.  However,  the  three  appen¬ 
dices,  on  mesh  relations  for  sequences  with  two-level  autocorrelation, 
on  the  shift-and-add  relations,  and  on  filter  integration  of  m-sequencer, 
present  material  that  to  the  author's  knowledge  is  original 

1.  INTRODUCTION 

During  the  past  decade,  a  new  class  of  modulation  functions  has 
been  developed  and  applied  to  communication,  guidance,  and  radar  sys¬ 
tems.  Called  pseudc-nolse  or  pseudo- random,  these  functions  are  dig¬ 
ital  in  nature.  7n  the  application  of  these  functions,  a  transmitted 
signal  having  a  small  number  of  modulation  states,  usually  two,  is 
employed.  The  signal  is  then  alternated  between  the  states  as  pre¬ 
scribed  by  the  digital  code. 

Compared  with  conventional  pulse  or  CW  radar,  radar  systems  using 
these  functions  may  exhibit  new,  and,  in  many  ways,  better  performance 
characteristics  such  as 

(a)  large  average-  to  peak-power  ratio 

(b)  unambiguous  measurement  of  range  to  large  ranges 

(c)  unambiguous  measurement  of  velocity  (Doppler  frequency!  to 

high  velocities 

(d)  fine  range  resolution 

(e)  fine  velocity  resolution 

(f)  resistance  to  both  sophisticated  and  power  jamming 

(g)  a  signal  difficult  ior  an  enemy  to  detect  because  of  its 

peak-free,  noise-like  spectrum. 

Ordinary  CW  radar  is  completely  satisfactory  with  respect  to 
average  power  and  measurement  of  velocity  but  neither  resolves  nor 
measures  target  range  satisfactorily.  On  the  other  hand,  a  pulse 
radar  using  short  pulses  and  a  low  duty  cycle  eliminates  the  rang¬ 
ing  problem  but  sacrifices  average  power  and  the  unambiguous 


measurement  of  velocity.  The  velocity  measurement  becomes  ambiguous 
because  many  short  pulses  must  be  used,  and  Doppler  must  be  deter¬ 
mined  through  a  pulse-to-pulse  phase-shift  measurement  instead  of 
through  an  intrapulse  frequency  measurement.  In  the  pulse-to-pulse 
method,  the  phase-shift  measurement  is  unambiguous  only  from  -71  to 
+77;  hence,  the  magnitude  of  the  largest  unambiguous  Doppler  frequen¬ 
cy  is  one  half  of  the  pulse  repetition  frequency.  The  result  is 
that  only  low  frequencien  can  be  measured  unambiguously. 

That  digital  modulation  functions  can  mitigate  the  disadvantages 
of  pulse  radar  to  some  extent  may  be  illustrated  as  follows.  Assume 
that  the  waveform  nay  be  transmitted  in  one  of  two  modulation  states, 
Ml  or  M2,  which  may  be  phase  states,  frequency  states,  amplitude 
states,  or,  in  the  limit,  off-on  states  (M2  =  MO).  For  a  system 
using  phase  modulation,  the  more  usual  case,  the  waveform  in  state 
Ml  is  illustrated  in  figure  1(a);  in  states  Ml  and  M2,  in  figure  1(b). 
In  this  case,  Ml  and  M2  differ  by  180  deg  in  phase. 

The  unambiguous  range  of  the  pulse  radar  may  now  be  increased 
by  transmitting  the  successive  pulses  in  either  state  Ml  or  M2  as 
prescribed  by  a  binary  sequence,  either  specified  or  random.  If  the 
sequence  is  periodic  with  a  period  of  L  digits,  the  unambiguous  range 
has  been  extended  by  a  factor  of  L.  Unfortunately,  the  integration 
time  has  been  extended  similarly.  If  the  sequence  is  aperiodic,  a 
more  complicated  situation  arises,  to  be  treated  later. 

If,  instead  of  desiring  to  increase  the  unambiguous  range,  we 
desire  to  increase  the  range  resolution,  and,  if  the  change  of  state 
can  be  accomplished  sufficiently  rapidly,  then  we  may  code  each  pulse 
with  a  binary  sequence  (fig.  1(c))  .  In  this  case,  the  range  resolu¬ 
tion  and  bandwidth  are  both  increased  by  a  factor  L. 

Alternately,  if  the  average  power  is  to  be  increased,  pulses 
modulated  as  in  figure  1(b)  may  be  concatenated  as  in  figure  1(d). 

The  limit  of  such  a  concatenation  is  modulated  Ci,'.  In  addition,  it 
can  be  shown  that  the  velocity  resolution  is  best,  in  a  certain  sense, 
if  angle-modulated  CW  is  used  (ref  16,  p.  24). 

If  immunity  to  countermeasure  is  prime,  then  an  aperiodic  sequence 
or  one  of  long  period  and  complicated  structure  might  be  employed. 

With  such  a  sequence,  it  will  be  difficult  for  a  jammer  to  produce  an 
"advanced"  waveform;  i.e.,  one  designed  to  deceive  the  receiver  by 
falsely  registering  as  a  ^er  target,  because  the  sequence  is  theoret¬ 
ically  or  practically  unpredictable.  In  this  case,  of  course,  the 
receiver  would  be  much  more  complicated, 

Ai ‘hough  it  would  be  desirable  to  hide  the  modulated  signal  in 
a  peak-free,  noise-like  spectrum,  it  is  usually  not  possible  to  do  so 
in  a  radar  system,  because  the  transmitted  power  is  so  large  as  to  be 
easily  recognized  even  with  substantial  spectrum  spreading.  The 


technique  may  be  useful,  however,  in  a  secret  communications  system 
where  power  may  be  minimized.  Modulation  of  the  message  by  a  binary 
sequence  i'  then  a  form  of  scrambling  that  spreads  the  spectrum 
(fig.  2(a)).  If  the  receiver  can  remove  the  binary  modulation 
through  knowledge  of  the  sequence,  the  received  spectrum  (message) 
is  collapsed  (fig.  2(b)).  At  the  same  time,  any  CW  present,  which 
may  be  CW  jamming,  is  spread,  while  uncorrelated  noise  is  not  col¬ 
lapsed.  The  message  may  then  be  separated  by  a  narrow-band  filter. 

Thin  illustration  of  the  application  of  various  modulation  func¬ 
tions  to  pulse  radar  has  thus  shown  how  certain  disadvantages  of  con¬ 
ventional  pulse  radar  may  be  alternatively  alleviated,  but  not  with¬ 
out  obtaining  certain  other  disadvantage::' .  The  illustration  implies, 
however,  that  by  proper  choice  of  the  modulation  function,  it  may  be 
possible  to  obtain  a  number  of  advantages  simultaneously  without  ob¬ 
taining  serious  disadvantages.  The  objective  of  this  report  is  to 
examine  the  various  modulation  functions  with  this  in  mind,  to  weight 
their  relative  advantages,  and  ideally,  to  arrive  at  a  best  modula¬ 
tion  function. 

2.  THE  RADAR  RECEIVER  AND  TRANSMITTED  WAVEFORM 


The  most  that  can  be  required  of  a  radar  receiver  is  that  it 
compute  a  probability  density  for  a  target  at  each  range  associated 
with  return  delays  r  and  each  velocity  associated  with  Doppler  fre¬ 
quencies  l^.1  Assuming  a  single  point-target  and  Gaussian  noise,  it 
has  been  shown  (ref  2,3,4)  that  no  receiver  can  extract  more  informa¬ 
tion  from  the  received  signal  than  a  correlation  (matched-filter) 
receiver.  For  a  transmitted  waveform  w(t),  the  delayed,  Doppler-shifted 
return  is 


r(t)  ss  w ( t  -  t)  +  n(t)  (1) 

where  n(t)  represents  additive  noise.  A  correlation  receiver  is  one 
that  computes  the  envelope 


env  C(  T,(J)  =  |T(t,4)  | 

=  ||  f  r*(t)  w(t  -  r)  e2771^  dt  |  (2) 

*1 

1  If  the  radar  system  becomes  ambiguous  in  t.  and/or  <^,  no  attempt 
is  made  to  resolve  these  ambiguities.  Therefore,  sufficiently 
strong  returns  from  the  ambiguous  (t,$)  region  will  be  falsely  inter¬ 
preted  as  returns  from  the  unambiguous  region  and  their  probabilities 
so  computed. 
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Figure  2(b).  Spectrum  after  code  demodulation. 
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where  I  is  the  interval  of  integration,  i.e.,  tl»e  conjugated 
signal  is  multiplied  by  a  delayed,  frequency-shifted  reference, 
integrated,  and  the  envelope  taken  (to  discard  rapid  fluctua¬ 
tions  in  C(t,<^)  that  carry  no  information).  Various  receivers 
of  this  type  have  been  designed1  and  are  described  in  the  refer¬ 
ences  . 

Although  a  correlation  receiver  extracts  all  available  in¬ 
formation  from  a  given  signal,  some  transmitted  signa’s  have 
more  suitable  characteristics  than  others .  The  function  that 
measures  the  suitability  of  a  waveform  in  this  respect  is  the 
signal  component  of  r(r,<^)  which,  after  a  change  of  origin,  is 
Woouward's  famous  ambiguity  function 

X(r,(^)  =  J  w(t)  w*(t  -  r)  e"2711^  dt  (3) 

I 

For  most  radar  purposes,  a  good  waveform  w(t)  is  one  for  which 
the  envelope  |X(T,<^)|  of  X(t,<^)  has  a  sharp  peak  at  the  origin 
and  is  small  elsewhere. 

We  expect  to  use  a  carrier  modulated  by  a  low-pass  modula¬ 
tion  function: 

w(t)  =  x ( t )  e2711^1  (4) 

Then , 

|X(t,^)  =  jj  x(i)  x*(t  -  r)  e  2vi^t  dt|  (5) 

and  so  the  pertinent  ambiguity  function  is  that  of  the  modulation 
waveform. 

Along  the  t  axis,  |x j  is  just  the  envelope  of  the  autocorre¬ 
lation  function  of  x(t): 

|X(T,0)|  =  |J  x(t)  x*(t  -  t)  dt  |  (S) 


1  This  kind  of  correlation  receiver  is  to  be  distinguished  from 
another  class  of  distance-measuring  correlation  systems  that  does 
not  delay  the  reference  signal,  and  that  essentially  tttces  the 
magnitude  of  env  at  the  origin  as  a  measure  of  the  nearness 

of  the  peak  of  env  C(t,$)  to  the.  ~rtjin.  Such  a  system  is  described 
by  B.M  Horio^,  "‘'oit-e-Modulated  Distance  Me.ju-ing  Systems,"  IRE 
Proceedings,  Vol  47,  1959,  pp  821-828. 
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For  certain  important  modulations,  the  ambiguity  function  is 
much  better  (i.e.,  smaller)  along  the  ~  axis  than  elsewhere  in 
the  plane.  The  modulation  period  (and  the  integration  interval) 
is  made  shorter  than  the  shortest  expected  Doppler  period  in  rder 
to  guarantee  that  correlation  of  the  received  signal  with  a  de¬ 
layed  but  not  f requencv-shif ted  signal  will  give  a  small  off  peak 
value. 

The  ambiguity  functions  oi  common  waveforms  such  as  pulsed 
sine  waves  and  "chirp"  signals  have  been  calculated  and  found  to 
be  far  from  optimal  (ref  3,  4,  30).  Except  for  a  few  modulations 
aa  yet  unexploited  in  systems  (ref  7),  the  study  of  better  modula¬ 
tion  functions  has  been  restricted  to  time-quantized  "telegraph 
signals"  like  those  in  figure  1.  Henceforth  we  will  consider  only 
functions  of  this  typer  the  modulator  chooses  among  a  fixed, 
finite  numbe.  ox  modulation  states  at  multiples  of  a  basic  time 
interval.  The  modulation  waveform  will  be  represented  by  a  se¬ 
quence  of  complex  numbers,  each  modulation  state  being  symbolized 
by  one  member  of  me  sequence. 

3 .  BINARY  DIGITAL  MODUMTION  OF  PULSE  RADAR 

To  begin  with,  we  will  investigate  the  simplest  case  of  binary 
modulation  of  pulse  radar,  where  there  are  two  modulation  states. 

It  is  sometimes  convenient  to  take  these  as 

a  =  +1  or  -1  (7) 

k 

ar.a  at  other  times  as 

ok  =  0  or  1  (8) 

Generally,  these  will  be  associated  in  the  order  shown  above;  i.e., 

a  =  1  -  2b^  (9) 

on  those  occasions  when  a  change  of  variable  is  performed.  If 
these  sequences  of  or  b^  are  associated  directly  with  the  modula¬ 
tion  waveform  via 

x ( t )  =  a  or  > ( c)  =  b  (10) 

K  K 

kt  <  t  ..  (k  +  1;. 
o  —  -  o 

then,  the  modulations  are,  respectively,  0  or  180-deg  digital  phase 
modulation,  or  on-off  digital  amplitude  modulation.  When  discus¬ 
sing  the  mathematical  properties  of  binary  sequences,  however,  the 
actual  form  of  modulation  will  be  disregarded  during  the  analvsi.  , 


and  the  above  change  of  variable  made  whenever  it  is  mathematically 
convenient.  Although  the  only  operation  between  a^'s  will  be  or¬ 
dinary  multiplication,  the  b^'s  will  be  either  multiplied  or  sub¬ 
jected  to  "addition  modulo  two,"  whichever  happens  to  be  convenient 
to  the  calculation  or  proof  at  hand.  The  latter  operation  is  defined 
so  that  it  corresponds  to  multiplication  of  the  ak's  and  is  symbol¬ 
ized  by  @  :  > 


x  +  -  (+)  0  3 

+  +  -  0  0  1 


-  -  *  110 
Multiplication  of  ±  1  Addition  modulo  two  of  0,1 


that  is, 


G  ©  1  =  1  =  1  ©  0 

(11) 

O'©  0  =  0=  1  ©  1 

(12) 

The  first  class  of  binary  codeu  waveforms  that  we  will  examine 

is  exemplified  by  figure  1(c).  In  these,  a  coded  putse  is  followed 

by  an  off-time  longer  than  the  pulse.  An  investigation  of  the  en¬ 
tire  ambiguity  function  will  not  be  made,  since  it  happens  that 

the  autocorrelation  function  itself  (which  is  the  zero-velocity 
cross  section  of  the  ambiguity  function)  is  not  one  particularly 
suitable  with  long  sequences. 

3.1  Coded  Words  and  Phase-Coded  Pulses 


If  the  two  modulation  states  Ml  and  ,V2  differ  only  in 
phase,  and  in  that  by  180  deg,  the  coded  pulse  Is  most  naturally 
represented  by  a  finite  sequence,  "word"  of  ak's.  The  autocorre¬ 
lation  function  of  a  coded  word  is  givt  ,  by 

L~t 

C  (t)  =  )  a.  a  ( 

a  L,  k  k  +  t 

k-.l 

where  t  is  a  discrete  variable,  r  -  *••»  -2,  -1,  0,  1,  2,-  •'< 

L  is  the  word  length,  and  the  autocorrelation  is  onnormalized ;  e.g., 
for  the  coded  word  (+1,  +1,  +1,  -1,  +1),  the  pulse  waveform  Is 


i 

» 

_ + 

1 
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and  aj  =  *2  =  *3  =  a5  =  +1 ;  a4  =  •  Then  Ca(l)  Is  co«puted  from 

the  partial  overlap 


+  +  +  -  + 

+  +  +  -  + 


and  is  equal  to 


C(l)=l+1-1-1=0 

a 


The  autocorrelation  function  for  this  pulse  is 


(14) 


and  furthermore,  this  is  representative  of  the  class  of  coded 
words  with  the  best  autocorrelation  function: 


C  (t)  = 
a 


j  0  or  + 

l" 


1 


T  =  ±  1,  ±  2, 

T  =  0 


±<L  -  1) 


(15) 


Unfortunately,  the  longest  known  such  "perfect"  word  is  only 
13  digits  long.  There  are  no  longer  ones  of  odd  length  (rof  33). 
The  existence  of  longer  perfect  words  of  even  length  hinges  on 
the  unsolved  problem  of  the  existence  of  long  "perfect"  periodic 
sequences,  which  will  be  mentioned  later.  Their  existence  is 
consldeied  unlikely.  Some  of  the  knowu  perfect  word?  are  (ref  1) 


L  Perfect  Word 

2  +  - 

3  +  +  - 

4  +++-,++-+ 

5  +  +  +  -  + 

7  +  +  +  --+- 

11  +  +  +  - 

There  has  been  some  experimental  investigation  oi 
words:  jc#<T)j  =  2  ,.r  3  off-peak. 
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"nearly  perfect" 


3.2  Coded  Vowels  and  On-Off  Coded  Pulses 


If,  on  the  other  hand,  the  two  states  Ml  and  M2(=M0) 
are  on  and  off,  the  coded  pulse  is  best  represented  by  a  finite 
sequence,  or  "vowel"  of  b^'s.  Such  pulses  have  been  used  in 
the  Venus  ranging  experiment  (ref  1,  16,  55).  On-off  coding  was 
used  because  the  return  was  not  phase-coherent. 

The  autocorrelation  function  for  a  coded  vowel  is 
L-t 

C  ( T )  =  V  b  b 
b  k  k+r 

k=l 


i.e.,  an  overlap  product  analogous  to  word  correlation.  The  per¬ 
fect  autocorrelation  in  this  case  is  defined  as 


C(t) 


f  N±i 

I  2  > 

i  0  or  1, 


An  example  of  a  perfect  vowel  is  11001011 ,  few  are  known. 
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4.  BINARY  CODES  FOR  CW 


The  case  of  continuous  transmission  remains;  we  now  treat 
digital  coded  CW  (the  conventional  but  somewhat  anomalous  term) 
and  in  the  next  section  restrict  the  codes  to  binary  sequences. 
Most  of  the  systems  built  so  far  employ  0-  to  180-  deg  phase 
modulation  (PM),  although  a  few  use  AM  or  FM.  PM  will  be  assumed 
in  most  cases,  and  accordingly,  the  sequences  will  correspond 
directly  to  the  modulation  when  written  in  terms  of  the  Most 

sequences  considered  will  be  periodic,  but  first,  an  aperiodic 
sequence  with  certain  advantages  will  be  treated. 

4.1  The  Fair  Coin  Sequence 


If  the  sequence  is  chosen  by  flipping  a  fair  coin,  it 
has  the  advantage  of  being  unpredictable  to  the  jammer.  There  is 
no  possible  way  to  jam  it  intelligently,  and  recourse  must  be  made 
to  brute-force  techniques.  Since  the  sequence  is  aperiodic,  no 
range  ambiguities  are  built  in. 

A  typical  section  of  a  telegraphic  modulation  function 
obtained  bv  actual  coin-fli oping  appears  as 


1  It  is  to  t y?  noted  that  the  association  previously  given  n  eq  (9) 
(  +  1  -*•  0)  and  (-1  -*•  +1)  does  not  carry  perfect  vowels  into  perfect 
words . 
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HEADS 


X(t) 

I  n 


TA'LS  1 


■4 - I - *~ 


i  I  I — h 


-►t 


The  normalized  autocorrelation  function  of  the  modulation  is 


lim 

C(t)  =  T  oo 


2T 


/ 


-T 


x(t)  x(t+t)  dt 


(18) 


The  autocorrelation  function  is  zero  except  for  t  =  0  because 

there  are  equal  numbers  of  agreements  and  disagreements  between 

a.  and  a.  .  The  former  contribute  +1  to  C(r) ;  and  the  latter, 

**  K41  T 

-1.  In  fact,  an  alternate  definition  of  the  normalized  autocor¬ 
relation  function  of  any  a^  sequence  is 


C  (r) 
a 


A-D 

A+D 


(19) 


where  A  is  the  number  of  agreements  during  the  period  of  corre¬ 
lation,  and  D,  the  number  of  disagreements.  In  the  limit  of 


equation  (18), 

then,  it 

is  clear  that 

/ 

i  1 

,  T  *  0 

C  (t) 
a 

1° 

>  1 T  |  >  1 

for  the  random 

sequence 

(20) 


C0(r) 


Now,  this  is  the  same  as  the  autocorrelation  of  a  single  rectan¬ 
gular  pulse.1  And  according  to  the  Wiener-Khi ntchine  Theorem,  the 
Fourier  transform  of  the  autocorrelation  function  is  just  the 
power  spectrum.  Hence,  the  power  spectrum  of  the  modulation  wave¬ 
form  x(t)  is  exactly  the  same  as  that  of  the  pulse,  which  is  the 

,  ...  sinsx 

familiar  (-— — ) . 


The  above  autocorrelation  function  is  not  practical  for 
radar  applications  in  that  by  its  definition,  it  requiree  an 

5  Although  *e  have  considered  Ca(t)  only  for  integral  ,  t  he  defin¬ 
ition  (  1 H  >  may  ‘«e  emploved  to  rtevelnn  the  i-onttnuous  curve  shown. 

IS 
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infinite  correlation  time  (ref  51).  If  the  correlator  integrates 
for  only  L  digits,  the  autocorrelation  is  represented  as  in  figure 
3;  for  a  given  t,  the  value  of  the  finite  sum  could  lie  any¬ 
where  between  the  dotted  lines  with  probability  1/2  it  will  lie 
within  the  shaded  region  if  L  is  large  (ref  6);  its  expectation  is 
the  heavy  solid  line,  i.e.,  the  autocorrelation  function  previously 
computed  for  the  limit  L  -*  cd  . 


Figure  3.  Finite  integrc cion-t ime  autocorrelation  for  the 
fair-coin  sequence. 

Lytle  (ref  32)  has  analyzed  the  more  restricted  case  of 
randomly  chosen  periodic  sequences. 

4.2  Pseudo-Random  Sequences 


There  exist  preassigned  (and  in  fact  periodic)  sequences 
having  autocorrelat *  on  properties  nearly  as  good  as  those  of  random 
sequence*  without  the  disadvantage  of  random  fluctuations  in  their 
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correlation.  They  exhibit  properties  close  to  those  that  random 
sequences  have  only  on  the  average.  For  this  reason,  they  are 
called  "super-random"  or,  more  commonly,  "pseudo-random."1 

The  following  postulates  obtaining  pseudo-randomnesu 
have  teen  often  used  in  various  combinations: 

PI)  In  each  period,  the  number  of  +l's  is  nearly  equal 
to  the  number  of  -l's. 

P2)  In  each  period,  half  of  the  runs  (groups  of  consec¬ 
utive  digits  of  the  same  kind)  are  of  length  one, 
one  quarter  are  of  length  two,  one  eighth  of 
length  three,* •••2  of  length  m;  etc.  For  a 
finite  period,  this  must  end  somewhere,  viz.,  at 
least  at  the  time  2-m  =  l/o. 

P3  For  a  period  of  length  2r,  every  n-tunxe  appears 
once  per  period. 

P4)  The  ambiguity  function  of  *he  corresponding  modula¬ 
tion  should  have  a  sharp  peak  at  the  origin,  *  ad  be 
small  and  uniform  to  some  specified  degree  else¬ 
where.  In  particular,  the  autocorrelation  shall  be 
two  level: 


VT)  - 1 


a  a 
k  k+T 


I  L,  t  =  0 

1  M,  T  4  0 


This  autocorrelation  is  cyclic,  unlike  the  pulse  case  previously  con¬ 
sidered;  e.g.,  for  the  sequence  (---++-+),  C  (2 )  is  computed 
from  the  comparison  of  the  sequence  with  its  cyclically  shifted 
■elf: 


+  +  -  + 


---  +  +  -  + 


where  |  represents  an  arbitrary  division  into  periods.  Here, 

C  (2)  =  +1  -1  -1  -1  +1  +1  -1 
a 


1  Generally  speaking,  a  sequence  is  called  "random"  if  the  method 
of  generation  has  no  known  properties  (e.g.,  natural  laws  or  pre¬ 
ferred  initial  conditions)  that  would  cause  departures  from  -an- 
domness.  a  sequence  is  called  "pseudo-random”  if,  regardless  1 
Its  raethoa  of  generation,  it  satisfies  some  particular  set  of  cri¬ 
teria  derived  from  the  expectation  values  of  a  random  sequence. 

As  a  result,  a  pseudo- random  sequence  usually  appears  more  random 
than  almost  all  random  sequences  of  the  same  length. 
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and  it  can  be  verified  that  this  particular  sequence  satisfies  P4 
but  not  P3: 


C  (T)  = 
a 


7,  r  =  0  (mod  7) 
•1,  r  4  0  (mod  7) 


(23) 


Sequences  satisfying  P3  are  called  de  Bruijn  sequences 
(ref  9,  10,  13,  23),  They  satisfy  PI  and  P2  but  do  not  satisfy  P4 
if  n  >  1,  as  will  be  shown  later.  Adding  one  +  to  the  above  exam¬ 
ple  forms  a  de  Bruijn  sequence  of  degree  n  =  3: 


-  +  +  +  - 


P3  is  now  satisfied  but  P4  is  not: 


C  (3)  =  -4 
a 


A  de  Bruijn  sequence  of  arbitrary  degree  n  can  be  formed 
by  the  following  rule:  let  (a,  to  a)  be  (-)  and  continue  the  se¬ 
quence  writing  a  (+)  whenever 
ing  ea 
yields 


It 


does  not  complete  an  n-tuple  appear¬ 
ing  earlier  in  the  sequence  (ref  9).  Thus  for  n  =  4,  this  algorithm 


—  +  +  +  + 


+  +  --  +  -  + 


It  is  not  possible  in  general  to  form  a  sequence  satisfying  P4  from 
a  de  Bruijn  sequence  by  omitting  the  last  element. 

4 . 3  Sequences  with  Two-Level  Autocorrelation 

Only  ?4  s  of  direct  interest  for  modulation  sequences. 
The  requirement  that  Ca(T)  be  two-level  has  certain  consequences,  a 
number  of  which  will  be  stated  1  .i  this  section  in  order  to  give  the 
reader  a  familiarity  with  this  class  cr  sequence. 

Consequence  I — First  of  all,  it  is  obvious  that  for  any 
binary  sequence,  the  unnor.malized  autocorrelation 


C  (t) 
a 


k=l 


a  a 

k  k+T 


(24) 


(or  in  iact  any  form  of  complete  pairing,  cyclic  or  not)  consists 
of  L  terms  each  of  which  is  (-t-1)  or  (-1),  and  hence, 


C  (t) 
a 


L  (mod  2)  for  all  t. 


(25) 


But  the  stronger  statement  is  true  for  autocorrelation  or  any  com¬ 
plete  pairing  of  binary  sequences  that 
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C  (t)  =  L  (mod  4)  for  all  t 
a 


<2G) 


This  can  be  seen  from  a  simple  example: 

If  any  group  of  Republicans  and  Democrats  is  paired 
off  in  any  fashion  against  another  group  of  the  same  size  and 
political  constitution,  an  even  number  of  political  arguments 
will  ensue,  equally  divided  between  Republicans  versus  Democrats 
and  Democrats  versus  Republicans. 


Proof:  Let  the  number  of  ordered  pairs  of  each  of  the 

four  kinds  be  symbolized  in  an  obvious  way:  N  is  the  number  of 

DR 

Democrat-Versus-Republican  pairings,  etc.  The  number  of  Demo¬ 
crats  in  each  group  is 


DD 


+  N_ 


DR 


DD 


+  N 


RD 


Hence,  the  number  of  disagreements  is 


(27) 


D  =  N  +  N _ =  2N. 


DR 


RD 


DR 


Q.E.D. 


(28) 


and  so  the  correlation  of  political  views  is  congruent  (mod  4)  to 
the  size  of  either  group: 


C  =  A  -  D 
a 


=  L 


4N. 


DR 


(29) 


Thus , 

as  stated. 


C  ( t )  =  L  ( mod  4 ) 
a 


(30) 


Consequence  II — Denote  by  N+  the  number  of  (  +  l)'s  in 
each  period,  and  similarly,  for  N_ .  Then, 

L 

)  C  (T)  =  L  +  M(L-l) 

L  a 


T=1 


L  L 

=  \  y  a  a 

L  L  k  k+ 

t=1  k=l 


r  L  -^2 

=  \l  ak 

Lk=i  J 


(31) 


L  +  M(L-l)  =  (N  -  S  ) 
+ 


2 


(ref  23) 


(32) 


Combining  this  with 


N  +  N  =  L 
+ 


(33) 


we  obtain 


N+  =  (L  ±  sTl  +  M(L-l)/2 


(34) 


where  the  signs  may  be  taken  in  either  order.  M  is  invariant 

under  a  ->  (-a  ).  Also,  since  (L  +  M(L-l))  must  bo  a  square, 
K  K 


and  hence,  for  L  >  2 


(35) 


M  >  -  1 


(36) 


For  M  = 


1, 


N  =  N  ±  1 
+ 


(37) 


These  are  the  binary  sequences  most  widely  used  for  modulation. 

The  fact  that  this  disparity  cf one  between  N  and  N  is  the  smallest 
attainable  for  a  two-level  sequence  with  L  %  2  follows  from  equa¬ 
tion  (34).  As  p.  result,  no  balanced  sequence  (N+  =  N_)  with  period 
greater  than  two  can  be  two-level;  in  particular,  no  de  Bruijn  se¬ 
quence  except  (+-)  is  two  level. 

Consequence  III — If  a  two-level  sequence  is  sampledevery 
q  digits,  wnere  q  is  prime  to  L,  the  sequence  so  produced  is  also 
two-level  with  the  same  levels  L  and  M  (ref  15).  This  can  be  seen 
by  sampling  the  autocorrelation  summation  of  the  original  sequence 
at  the  same  rate,  for  each  value  of  t. 

Consequence  IV- -The  corresponding  sequence  of  (b^  =  C  or  1) 
also  has  two-level  autocorrelation: 


C  (t) 

a 


=  Y  (1  -  2b)  (1  -  2b  ) 

L  k  k+T 

k=l 

=  4C  (r)  +  N+  -  3N_ 


(38) 


CW(T>= 
b 


!(m  +  ; 


T  =  0 


(39) 


3N  -  N  )/4  ,  r  ±  0 
+ 


Each  run  of  r  "ones”  contributes  r  to  C  (0);  it  contributes  (r-1)  to 

C  (1).  Hence  the  number  of  runs  of  (l)’s  in  each  period  is 
b 
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c-!0) 


c,  (O)  -  C,  (1)  =  N  -  (M  +  3N  -  N  )/4 

»J  13  —  — 

=  (L  -  M)/4 


and  this  must  also  be  the  number  of  runs  of  (0)'s,  since  the  two 
kinds  alternate.1  Each  run  of  (r  >  1)  "ones"  contributes  (r  -  2) 
to  C^(2),  and  each  run  of  a  single  "one"  contributes  zero,  with 
the  exception  that  at  each  run  of  a  single  zero,  there  is  an  extra 
contribution  of  one.  If  the  number  of  runs  of  a  single  zero  is  de¬ 
noted  by  and  of  a  single  one  by  JB^, 


C,  (0)  -  Cu (2)  =  2((L  -  M)/4  - 
b  b 


B1>  +  1 


(41) 


°B  +  1B1  =  (L  -  M)/4  (42) 

i.e.,  half  of  the  total  number  of  runs  of  symbols  of  either  kind  are 
of  length  one.  Hence  P4  guarantees  the  first  clause  of  P2.  No  such 
simple  statement  is  possible  about  the  other  clauses  of  P2 . 

Consequence  V — It  is  sometimes  necessary  to  know  the  number 
of  ordered  pairs  of  each  kind  of  (a  a  )  appearing  in  the  C  (t)  sum¬ 
mation.  Denote  these  bv  N  ,  N  ,  N  .and  N  ,  and  let 

++  +- '  -+’  — * 

m  =  (M  +  3N_  -  N+)/4  (43) 

Applying  the  reasoning  of  the  political  example  to  Cb, 

N _ =  m 

-  m  (44) 

-  m)  . 


Consequence  VI — The  autocorrelation  function  of  a  two- 
level  sequence  is  the  same  as  that  of  a  periodic  rectangular  pulse 
train  except  for  a  possible  change  in  d-c  level.  By  the  Wiener- 
Khintchine  Theorem,  the  power  spectrum  must  therefore  be  the  same 
also.  It  is  a  line  spectrum  with  onvelope  (sin8x/xs),  and  there 
is  a  d-c  value  appropriate  to  the  particular  value  of  M. 

If  a  sequence  is  used  to  phase-modulate  a  carrier,  the 
spectrum  will  be  the  convolution  of  the  sequence  spectrum  and  the 

1  The  average  run  length  2L/(L  -  M)  Is  equal  to  two  (the  value  for 
a  random  sequence)  only  if  M  is  zero. 


s  =  N  =  N 
+-  -+ 


N  =  N  -  (N 
++  + 
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carrier  spectrum.  This  will  be  free  of  peak6  (and  thereby  facil¬ 
itate  hiding)  provided  the  sequence  has  a  small  d-c  level.  Hence 
for  the  purpose  of  hiding,  the  choice  of  M  =  -1  is  optimum,  since 
it  yields  the  most  nearly  balanced  sequence. 

In  most  systems  built  to  date,  the  sequence  period  is 
kept  less  than  the  ieast  expected  Doppler  period,  in  order  that 
only  the  region  of  the  (r,^)  plane  near  the  t  axis  is  accessible 
to  signals.  In  such  a  case,  the  optimum  choice  of  M  to  keep  [xj 
3mall  off-peak  wo  id  be  M  =  0.  Tie  re  is  exactly  one  such  sequence 
known  (-*++-)  (its  inversions  and/or  rotations  are  not  considered  to 
be  Independent).  There  is  no  known  proof  that  longer  sequences  with 
M  x  0  do  not  exist,  but  some  restrictions  can  be  put  on  the  possible 
periods.  From  equations  <32,  33),  it  follows  that  the  period  L  must 
be  an  even  square: 


L  =  4na  n  =  1,  2,  3,  •  •  •  (45) 

Turyn  has  shown  further  that  n  may  not  be  a  power  of  a  single  prime: 
n  ^  P01,  P  prime.  The  smallest  values  of  L  remaining  are  144  and  400 
(n  =  6  and  10,  respectively j . 

An  additional  relation  for  two-level  sequences  is  derived 
in  Appendix  A. 

4 . 4  Shift-Register  Generators 

4.4.1  General  I/x; l c.  (ref  23) 

A  shift-register  generator  (SRG),  of  degree  n  con¬ 
sists  of  n  fltp-flopu,  a  clock  pulse-generator,  and  a  logic  circuit, 
as  shown  in  figure  4.  For  each  flip-flop,  "off"  is  symbolized  by 


'OUTPUT 


Figure  4.  The  general  SHG. 
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(0)  and  "on”  by  (1).  When  the  clock  pulses  arri*  ,  the  contents 
of  each  flip-flop  are  transferred  to  the  next  st<ige,  as  shown. 

The  new  state  of  the  first  stage  is  an  ..’•bitrary  Boolean  function 
of  the  previous  contents  of  the  register.  The  function  specifies 
a  (0)  or  (1)  as  desired  for  each  ot  the  2n  states  of  the  SRG; 
hence,  there  are 


different  such 


fund  ions  . 


!f  the  logic  is  correctly  chosen,  the  SRG  may 
produce  a  long  sequence.  The  longest  possible  SRG  sequences,  for 
a  given  n,  are  the  de  Bruijn  sequences  (1,  -  2n)  ;  it  has  been  shown 
that  there  are 


different  de  Bruijn  sequences  of  degree  n.  As  previously  noted, 

none  of  them  are  two-level  if  n  >1. 

/ 

The  logic  uniquely  specifics  the  successor  to 
every  state  of  an  SRG.  But  although  each  state  has  a  unique  suc¬ 
cessor,  it  is  possible  for  a  Si.ate  to  have  cither  zero,  one,  or  two 
predecessors.  For  example,  consider  a  d -stage  SRG  for  which  the 

new  state  oi  the  first  stage  is  the  ordinary  arithmetic  product  of 
the  previous  states  of  the  last  two  stages: 


The  sequences  oj  states  tor  tins  register  an* 


Here,  (110)  and  (100)  have  no  predecessors;  (010)  and  (000)  each 
have  two  predecessors.  We  will  return  to  the  general  SRG  later. 


4.4.2  Linear  Logic  and  M-Sequencea (ref  12,  14,  15,  17,  18) 

fin  important  subclass  of  SRG's  consists  of  those 
linear  logic.  A  linear  logic  unit  is  one  which  computes  the 
sum  mod  2  of  the  contents  of  whichever  stages  feed  the  unit.  It  can 
easily  he  shown  that  every  state  of  a  linear  SRG  has  a  unique  prede¬ 
cessor,  and,  as  a  result,  the  SRG  decomposes  the  space  of  2  states 
into  a  number  of  closed  cycles.  Since  the  sum  mod  2  of  any  number 
of  zeroes  is  zero,  the  'all-zeroes"  state  always  forms  a  separate 
cycle.  If  the  stages  that  are  to  be  tapped  are  chosen  correctly, 
there  is  only  one  other  cycle,  of  period  (2n-l). 

Consider  a  (3,5)  linear  SRG,  which  is  one  such 
"linear-maximal"  (or  "linmax")  SRG: 


If  the  initial  contents  are  00001,  the  contents  after  successive 
clock  pulses  will  be 

(init'al  contents)  1)  00001 

2)  10000 

3)  01000 

4)  00100 

5)  10010 

6)  01001 

7)  10100 

8)  11010 

9)  01101 

10)  00110 
etc . 

28)  10101 

29)  01010 

30)  00101 

31)  00010 
1  =  32)  00001 

The  sequence  appearing  at  the  output  is 

. . .1000010010110011111000110111010. . . 

If  one  (0)  were  added  to  the  run  of  tour  (0)'s  the  sequence  would 
be  a  de  Bruijn  sequence.  Hence,  this  SRG  sequence  almost  satisfies 
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PI,  P2,  and  P3  exactly,  as  it  stands 
a  two-level  autocorrelation: 


But  most  important,  i 


31  =  L  ,  j  T  =  0 

Ca(T>  -  1  <46) 

-i  ,  [r  4  f' 

It.  is  not  difficult  to  prove  that  every  linear- 
maximal  sequence  (called  m-sequence)  has  a  similar  autocorrelation. 

In  the  first  place,  a  linear  logic  satisfies  superposition:  if 
two  identical  linear  SRG's  A  and  B  have  arbitrary  contents,  and  a 
third  identical  SRG  C  contains  the  sum  (mod  2)  of  the  contents  of 
A  and  B,  then  the  output  digit  from  C  will  be  the  sum  (mod  2)  of 
the  outputs  cf  A  and  B.  From  this  there  follows  an  important 


Shi ft -and -add  property — If  an  m-sequence  is  shifted 
by  r  digits  (r  neither  zero  nor  a  multiple  of  L) 
and  added  to  the  unshiited  sequence  mod  2,  the  sum 
sequence  is  the  same  sequence  shifted  by  some  other 
number  of  digits  r',  e.g.,  for  the  (3,5)  sequence 
and  t  =  2, 

, . . 100001001011001I.U1000110111G1010. . . 

. . .1000010010110011111000110111010. . . 

.  .  .lOOlOilOOlllllQOOUOl  1101010000 

t'  -  28 


The  dependence  of  r'  on  r  is  erratic  as  may  be  seen  from  the  sam¬ 
ple  sequence: 


T 

:  1 

2 

3 

4 

5 

6 

7 

8 

S 

10 

11 

12 

13 

14 

15 

t' 

:  14 

28 

5 

25 

3 

10 

16 

19 

24 

6 

23 

20 

30 

1 

22 

16 

17 

18 

IS 

20 

21 

22 

23 

21 

23 

26 

27 

28 

29 

30 

7 

18 

17 

8 

1.2 

27 

1.5 

11 

g 

4 

29 

21 

2 

26 

13 

This  dependence  is  not  given  completely  by  any  known  function;  its 
known  properties  are  stated  in  Appendix  B. 

Since  a  (O';  occurs  in  the  sum  sequence  above  if  and 
only  if  the  digits  in  the  addend  sequences  agree,  and  since  there  is  one 
more  (1)  rhan  (0)  In  a  period  of  an  m-sequence,  the  autocorrelation 
function  must  be 


C 

a 


(T> 


-  =  0 


(47) 


as  we  set  out  to  prove. 


Nothing  a. -ove  guarantees  the  existence  of  m-sequences, 
but  it  has  been  shown  (ref  15,  18)  that  there  are  in  fact  (<)(2!'‘-l))/n 
m-sequences  of  degree  n.  (<J>  is  Euler's  function:  (^(k)  is  the  number 
of  integers  less  than  k  which  are  prime  to  k,  including  1.)  To  find 
the  logic  (i.e,,  which  stages  are  tapped)  cf  these  linmax  SRG's  is  to 
find  the  primitive  irreducible  polynomials  over  the  Galois  field  of 
two  elements.  This  computational  problem  has  been  solved  with  a  dig¬ 
ital  computer  for  all  linmax  SRG's  with  n  not  greater  than  16,  and 
for  some  SRG's  of  each  n  up  to  34,  and  the  results  are  listed  in 
Appendix  C  of  reference  20.  Among  these  permissible  tap  combinations 
are  (1,2),  (2,3),  (3,4),  (3,5),  (5,6),  (4,7),  (4, 5, 6, 8),  (5,9), 

(7,10),  (9,11),  (6,8,11,12),  ...,  (20,33),  (7,32,33 ,34).  A  few 
longer  linmax  SRG's  are  known,  e.g., 

(126,127),  l,7,15,8091),...(2127-2,  2127-]). 

Given  one  linmax  SRG  of  degree  n,  all  others  of 
degree  n  may  be  found  from  the  Prime  Sampling  Theorem  previously 
stated,  which  in  this  case  at  least,  turns  out  to  be  exhaustive.  The 
output  is  sampled  every  qth  digit,  where  q  is  odd  and  prime  to  L.  As 
soon  as  (2n-l)  digits  have  been  so  produced,  the  SRG  that  will  pro¬ 
duce  this  new  sequence  can  be  constructed  by  solving  n  linear  equa¬ 
tions  for  the  tap  positions.  Repeating  the  process  produces  all 
linmax  SRG's  of  degree  n  (ref  17). 

Another  property  of  m-sequences  that  has  been  put 
to  use  in  FM  systems  (ref  25,  fig.  5)  is  that 


k-1 


>  (or  (1  +  b,  )  depending  on  the  choice  of  origin  k^l) 

i+s  i+s  (4g) 


where  the  summation  is  mod  2,  and  where  s  is  fixed  for  a  giver  sequence. 
It  can  easily  be  shown  that  this  property  is  true,  and  that  s  is  equal 
to  the  number  of  digits  separating  the  beginning  of  the  run  of  n  ones 
from  the  beginning  of  the  run  of  (n-1)  zeroes. 


The  sequence  of  a  linmax  SRG  can  he  obtained  with  any 
desired  delay  by  adding  (mod  2)  the  contents  of  some  combination  of 
the  n  stages  (ref  35).  This  could  obviously  be  very  useful  in  a  radar 
system,  in  order  to  produce  reference  signals  with  different  delays. 

4,5  Other  Sequences  with  M  =  -1 


There  are  two  other  known  classes  of  sequences  with  M  -  -1, 
known  as  "Perron"  (or  "Legendre"  or  'quadratic-residue")  sequences 
(ref  15,  26,  36)  and  "twin-prime"  sequences  (ref  26,  39),  rv  pec  ti  veiy . 

Perron  sequences  are  of  period  (4ra-l)  when  this  quantity  is 
a  prime.  Their  construction  is  illustrated  by  the  case  for  m  -  3 . 

The  residues  of  the  squares  of  successive  integers  modulo  (4(3)-l  =  ll) 
are 


0,  1,  4,  9,  16 (mod  11)  =  5,  25(mod  11)  =  3,  3,  5,  9,  4,  1,  0,  1,  4,... 
and  the  sequence  Is  defined  by 


+1  if  k  is  among  0,  1,  3,  4,  5  9 

ak=<[  '  (49) 

-  1  otherwise , 

i.e., 

+  +  -  +  +  +  --  --*•- 

Except  for  m  =  1  or  2,  these  sequences  are  not  m-sequences  and  do 
not  have  the  shif t-and-add  property.  This  is  an  advantage  in  cer¬ 
tain  respects,  as  will  be  seen  later. 

Any  periodic  sequence  can  be  generated  by  some  SRG,  and 
in  a  physical  system  this  is  often  the  method  most  economical  of 
components.  The  above  sequence  could  be  generated  by  any  one  of  a 
number  of  non-linear  SRG's  of  n  =  5  stages.  (Four  would  be  too 
few  because  +  ~  +  +  occurs  twice.)  However,  the  sequence  does 
not  contain  21  of  the  32  possible  quintuplets  of  (  +  ) 's  and  (-)’s. 

If  an  SRG  were  used  which  decomposed  the  set  of  all  quintuplets 
into  several  closed  cycles  (as  any  linear  SRG  must,  for  example), 
it  would  have  to  be  started  with  one  of  the  eleven  legitimate 
quintuplets,  and  if  an  error  occurred  might  jump  into  an  incorrect 
cycle  composed  of  some  of  the  21  "bad"  quintuplets .  As  previously 
noted,  non-linear  SRG's  need  not  generate  closed  cycles.  This 
fact  can  be  used  to  generate;  e.g.,  Perron  sequences  with  a  non¬ 
linear  SRG  which  avoids  the  above  two  difficulties,  with  each 
illegitimate  state  (here,  one  of  the  21  bad  quintuplets)  leading 
back  into  the  cycle  of  legitimate  states. 

The  twin-prime  sequences  have  period  L  =  p(p+2)  where  p 
and  (p+2)  are  both  prime,  and  are  formed  by  a  method  similar  to 
that  for  Perron  sequences.  The  twin-prime  sequence  for  p  -  3  is  the 
same  as  the  (1,4)  m-sequence.  The  p=5  sequence  in  notation  is 

loiiioooimioiiiooiooooioioncoioo 

The  remarks  on  the  generation  of  Perron  sequences  apply  also  to 
twin-prime  sequences,  of  course. 

•i.6  Acquirable  Codes 

To  locate  a  given  target  in  range  by  testing  every  inte¬ 
gral  value  of  t  wouid  require  L  tests.  However,  only  log2  b  tests 
would  be  needed  if  the  target  could  be  localized  to  within  1/2  L, 
then  to  within  1/4  L,  etc.  The  best  practical  solution  to  date  is 
to  use  certain  sequences  obtained  by  combining  shorter  sequences 
with  relatively  prime  periods,  so  that 


(50) 


P1  P2 


where  is  tne  period  of  the  i  ^  subsequence.  For  these  code9, 
at  most 


N  =  pi  +  p2  +  •  •  ■  +  Pk  (51) 

tests  are  needed,  and  on  the  average,  only  N/2.  These  codes  have 
been  used  in  the  Goldstone  deep-space  ranging  system  (ref  1,  49). 

4.7  The  Ambiguity  Function  for  Sequences 

If  the  period  of  the  modulation  waveform  is  not  kept 
small  in  comparison  with  a  poppler  cycle,  the  output  of  a  corre¬ 
lator  is  given  not  by  the  autocorrelation  function  but  by  the 
full  ambiguity  function. 

The  squared  magnitude  of  the  ambiguity  function  for  se¬ 
quences  is  defined  analogously  to  that  for  continuous  functions: 
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(52) 


where 


u  =  exp  (2iri/L) 


(53) 


This  is  obviously  equal  to  L2  at  the  origin  and  is  zero  elsewhere 

along  the  Doppler  frequency  (s)  axis.  Lerner  (ref  21)  has  computed 

C  for  m -sequences  and  found 
TS 


C 


Ts 


=  1 


,  t,s  =  0,0  (mod  L) 

,  T  =0  (mod  L) ,  s  4  0  (mod  L)  (54) 

,  t  ^  0  (mod  L) ,  s  =  0  (mod  L) 


=  L  +  1  ,  elsewhere 


i.e,,  a  peak  at  the  origin  and  a  plateau  elsewhere  cut  by  valleys 
along  the  axes.  This  calculation  uses  the  shi f t-and-add  property 
of  rr-sequences ;  other  (e.g.,  Perron)  sequences  with  two-level  auto¬ 
correlation  functions,  but  not  having  the  shi f t-and-add  property  do 
not  in  fact  exhibit  the  same  peak-free  plateaus. 


It  can  be  shown  that 


T,S 


(55) 


And  hence,  this  ambiguity  function  Is  very  nearly  the  smoothest 
obtainable,  i.e.,  has  the  smallest  maximum  value  o.  ..-peak;  the  only 
further  smoothing  possible  would  be  the  lowering  of  the  plateau 
height  from  <L  +  1)  to  L  by  sacrificing  the  valley  along  the  T_axls 


C  =  L 

YS 


=  L  ,  t,s  e  0,0 

=  0  j  elsewhere  on  the  a-axis  (necessary)  (56) 

o 

=  L  ,  at  the  other  (L  -  L)  points 


This  requires  that 


m  =  -nr 


and  from  equation  (A-26). 


(N  -  N  )  =  L  +  M(L  -  .1) 

+ 

*  n2  +  n(n2  -  1)  (5f 

i.e.,  (n^  +  -  n)  must  be  a  perfect  square.  This  is  never  the  case 

for  1  <  n  <  200,  at  least;  and  herce,  there  are  no  such  "smoothest- 
possible"  sequences  for  L  <  40,000. 

Since  moving  off  the  ;  axis  implies  a  degradation  of  the 
normalized  ambiguity  function  magnitude  from  1/L  to  *]  L  +  1  /  1/ 

n/TT~  for  an  m-sequence,  it  is  clearly  advantageous  to  keep  the  inte¬ 
gration  time  (Lt0)  much  less  than  one  Doppler  cycle.  If  the  integra¬ 
tor  is  a  filter  instead  of  some  type  of  block  integrator  (one  which 


computes  £  afc  a^  exactly)  however,  the  output  off-peak  may  be  de- 
1 

graded  in  two  ways.  First,  the  infinite  tail  of  the  filter  impulse 
response  will  remember"  some  input  contributions  extending  back 
over  past  Doppler  cycles.  By  using  a  short  time  constant,  this  may 
be  made  as  small  as  desired.  But  then  the  second  effect  appears: 
the  off-peak  autocorrelation  is  only  constant  if  the  integration  is 
a  block  integration  extending  over  an  integral  number  of  sequence 
periods.  Thus,  the  output  of  a  filter  will  fluctuate  appreciably 
unless  the  time  constant  is  much  longer  than  the  sequence  period. 
(This  effect  is  investigated  in  Appendix  C.)  Since  the  normalized 
off-peak  autocorrelation  (-1/L  for  an  m-sequence)  becomes  better 
with  increasing  sequence  period  L,  It  is  all  the  more  important 
that  the  bit  time  (tQ)  be  short. 

5.  BLOCK  DIAGRAM  OF  A  SYSTEM 


t  FI  shoe  in ,  and  Rittenbsch  have  described  several  system! 
that  use  m-sequ«nca  modulation  (ref  24,  25,  73).  Without  specify¬ 
ing  the  particular  kind  of  modulation,  consKier  the  following 
system: 
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Figure  5.  Block  diagram  of  a  system. 


Each  rotation  of  the  continuously  variable  phase-shifter 
causes  SRG  2  to  drop  one  digit  further  behind  SRC  1  in  the  m- 
sequence.  The  output  of  the  RF  demodulator  is  a  Doppler  signal 
chopped  by  the  m-sequence,  and  hence,  consists  of  components 
too  high  in  frequency  to  pass  any  of  the  filters,  which  serve 
as  the  integrators.  When  the  reference  code  from  SRG  2  has  the 
correct  delay,  the  code  demodulator  reassembles  the  chopped  Doppler 
into  an  unchopped  wave,  which  passes  some  filter.  Range  is  meas¬ 
ured  by  the  number  of  digits  SRG  2  has  to  drop  behind  SRG  1  to  pro¬ 
duce  a  peak  in  the  oul_,ut  of  some  filter. 

The  system  can  be  altered  in  certain  ways  to  utilize  trans¬ 
mitter  leakage  as  the  local  oscillator.  Unlike  the  above  system, 
this  alteration  requires  a  sequonce  with  the  shi f t-and-add  property. 
This  version  has  boon  Implemented  by  Craig  et  al  using  FM  and  5- 
stage  1 inmax  SRG's  (L  -  31).  They  report  "an  over-all  receiver 
sensitivity  of  -140  dbm  on  slowly ^moving  targets."  Since  the  sys¬ 
tem  must  search  serially  In  range,  the  acquisition  time  will  be  com¬ 
paratively  large. 

<i .  NON-BINARY  CODES  FOR  CW 

If  the  number  of  modulation  states  is  greater  than  two,  the 
multiplication  operation  and  the  autocorrelhtion  function  must  be 


:io 


$ 


redefined.  If  0,  120,  and  240  deg  PM  is  us  the  natural  representa¬ 
tion  of  the  consists  of  the  cube  roots  oi  unity,  and 


C  (t)  =  V  a*  a, 
a  l_  k  k+T 

k 


In  general,  q  states  of  equiangular  Pf.'  are  represented  by  the  qth  roots 
of  unity  and  this  autocorrelation  is  the  natural  one. 

p 

One  class  of  perfect  codes  is  known,  of  period  q  ,  for  any  q 
(ref  50).  One  period  of  one  of  these  codes  consists  of  the  following 
powers  of  s  (s  =  exp  (2Di/q)): 


1,2,3,  ...  q,  2,4,6,...,  2q,  3,6,9 -  3q , . ,q ,2q ,3q , . . . ,q‘ 


For  q  =  3  this  is  the  sequence 


2  2 

s,s  ,l,s  ,s, 1,1, 1,1 


The  autocorrelation  function  is 


C  <t> 

a 


L  =  q  ,  T  =  0 
0  ,  r  4  0 


The  limitation  on  length,  given  q,  is  a  serious  disadvantage  of 
these  sequences.  Sequences  of  any  desired  length  but  with  a  poorer  auto¬ 
correlation  can  be  obtained  from  multi-state  SRG's.  Here  q  may  be  any 
power  of  a  prime, 

m 

q  =  p  ,  (£ 


and  linear  SRG's  are  those  which  have  a  logic  consisting  of  multiplication 
by  0,1,2....,  or  (q-1)  at  each  tap,  followed  by  addition  mod  q.  As  before 
there  exist  liumax  SRG's  of  every  length  n  with  autocorrelation 


C  (t)  = 
a 


L  ■=  q  -  1  ,  T  =  0 
-i  ,  t  4  0 


Under  addition  mod  q,  these  sequences  have  the  shi f t-and-add  property. 
Connections  for  some  linmax  SRG's  are  given  bv  the  coefficients  of  the 
primitive  irreducible  polynomials  in  references  19  and  41;  e.g.,  for 

’  O 

q  -  3,  reference  19  gives  (x  -  x  -  2)  as  a  primitive  irreducible  poly¬ 
nomial.  The  coefficients  are  1, 0,-1, -2,  and  the  last  three  specify  a 
1  i  mnax  SRG : 
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The  m-sequence  produced  is  ...00101211201110020212210222... 

If,  on  the  other  hand,  the  modulation  states  are  "off,"  "on," 
0  deg  ,  and  "on"  180  deg,  the  natural  representation  is  a^  =  0, 

+1,  -1.  Tompkins  (ref  45)  has  found  all  of  the  perfect  codes 
of  this  kind  for  L  <  19  by  trial.  Fur  L  >  13,  these  have  few 
non-zero  terms,  a  distinct  drawback,  since  the  average  power 
suffers  accordingly. 

Non-binary  sequences  have  not  yet  been  used  in  radar  sys¬ 
tems  to  the  author's  knowledge. 
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APPENDIX  A.  MESH  RELATIONS  FOR  SEQUENCES  WITH  TWO-LEVEL  AUTOCORRELATION 


Let  (a  =  ±  1)  be  a  sequence  of  period  L  with  a  two  level  periodic 
autocorrelation.  If  q  ie  a  factor  of  L,  then  the  sequence  can  be  di¬ 
vided  into  q  "meshes"  by  sampling  every  qth  digit  beginning  with  any 
of  the  first  q  digits.  We  will  derive  relations  between  the  sums  Sq^ 
of  the  digits  in  the  various  such  meshes: 


qi  "  L  k 

k 


k  s  i  (mod  q) 


(A-l) 


One  tool  with  which  such  a  separation  into  meshes  can  be  accom-- 
plished  is  exp  (27liT/q): 

L  L-l 


^  C(r)  exp  (2-rriT/q)  =  L  •  1  +  M  ^  exp  (27liT/q) 


(A-2) 


T=1 

L  +  M(L  1),  q  =  1 
L  -  M  ,  q  >  1 


But  also, 

Li 

exp  (2i7iT/q)  =  ^  2.  ak  ak+t  SXP  (2l7iT/q* 


(A-3) 


t=1  k=l 


*k  exp  ( —  2tt i k/q )  )  t  exp  (27ii(k+T)/q) 


\  ak  exp  (-2771  k/q)  £  exp  (+2rir/q) 


a,  cos 
k 


(2tt k/q ) -  i  Y  sin(2irk)/q)  y  ak  cos(2nk/q)+i  ^  ak  sin(2t?k/q) 

1c  J  [_  k  k 


V  1 

2 

L  +  M(L  ‘  h  -1)  = 

)  cos(2rk/q) 

-  k 

)  sln(2T?k/q) 

where  6  is  the  Kronecker  delta. 

We  now  specialize  this  temporarily  to  the  case  M  =  0.  These  se¬ 
quences,  which  were  described  in  section  A. 3  have  periods  of  the  form 


38 


L  =  4n 


n  =  1,  2,  3,...  ,  n  |  p  ,  p  prime 

2 


Hence  L  is  always  divisible  at  least  by  2,  4,  n,  and  n 
For  q  =  2,  the  mesh  relation  above  yields 

41,2  -  (S22  -  V2 


(A-5) 


(A-6) 


We  know  also  that 

q 


Vs  .  = 
L  qi 


L 

r-1 


<  a,  =  N  -  N  =  2n 

L  k  + 


i=l  k=l 

Choosing  S22  >  S21  arbitrarily,  we  find  that 


(A-7) 


and 


o  —  n 
21  *  ° 


S22  «  2n 


(A-  8) 
(A-9) 


i.e.,  the  first  mesh  of  alternate  digits  is  balanced  and  the  second 

contains  the  entire  Imbalance  between  N  and  N 

+ 

For  q  =  4, 


4n  =  <S44  '  S42)  +  (S41  "  S43} 


From  the  case  q  =2,  we  have 


and 


S41  +  S43  =  ° 


S42  +  S44  =  2n' 


(A-10) 

(A-U) 
(A  12) 


Again  choosing  arbitrary  signs  and  indices  wherever  possible, 

1/2 


and 


"S<13  ‘  b41  ■  +  (S42  S44} 


S  „  +  S.  -  2n 
42  44 


(A- 13) 
( A-14  ) 


with  the  auxiliary  e.ndit.  ns  that 


‘  Tills  was  first  ai.  covered  and  proved  by  N.  Karav’anis  and  C.  A. 
Morrison  i.lcng  with  a  number  of  other  results  concerning  sequences 
with  M  -  0  that  are  not  included  here. 


(A-15) 


lS4i  I  5  n 


S42  >° 


S4i  >0 


S  . .  2  n(mod  2) 

4i 


(A-16) 

(A-17) 

( A— 1 8 ) 


For  any  n,  one  solution  is  given  by 


S43  ~  S41  “  S42  _  S44  “  n 


(A-19) 


for  even  n,  Lhere  is  also  the  solution 


S4!  =  S42  -  S43  '  °;  S«  =  2n  (II) 


(4-20) 


These  are  all  of  the  solutions  except  for  a  small  number  of  others 
that  may  be  found  straightforwardly  by  listing  the  squares  of  suc¬ 
cessive  composite  values  of  S  ,  factoring  these  into  S  and  S  , 
and  calculating  n  =  (S^2  +  S  )/2.  This  is  required  to  be  an  inte¬ 
ger  and  to  have  the  same  par_ty  as  S  ,  as  noted  above.  The  only 
such  irregular  solutions  for  n  <  20  are 


n  =r 

10; 

S4i  = 

(8, 

4, 

-8, 

16)  or  (6,  2,  -6,  18) 

n  = 

15: 

CO 

II 

3; 

-9, 

27) 

n  = 

20: 

S4i  " 

(12, 

4, 

-12, 

,  36) 

The  others  for  n  <  40  are  for  n  =  26(2),  30(2),  34(2),  35,  and  40(2) 

Relations  similar  to  the  above  can  be  written  for  each  of  the 
factors  of  any  sequence  with  two-level  autocorrelation.  The  ones 
above  ere  suited  to  a  computer  search  for  sequences  with  M  =  0. 

Such  a  search  might  begin  by  writing  the  admissible  combinations 
of  lengths  of  runs  of  (+l)'s  and  (-l)'s  from  the  two  rules  govern¬ 
ing  them;  e.g.,  for  n  =  2  these  are 

•  1  -fl 

f  5, 2, 2,1 
3, 1,1,1  J  ‘*.3,2,1 
l  3 , 3 , 3 , 1 


2. 2, 1,1  \ 


■to 


5 . 3 . 1 . 1 

4. 4. 1.1 


where  the  choice  has  been  made  that 


N  >  N 
+ 

These  would  then  be  interleaved  in  all  possible  independent  ways 
and  each  of  the  resulting  sequences  tested  cyclically  to  see  if  it 
obeys  the  mesh  relations  for  all  q  dividing  L.  Any  sequences  sur¬ 
viving  this  test  would  then  be  tested  for  two-level  autocorrelation, 
a  considerably  longer  process. 

We  now  return  to  tne  case  of  general  M.  For  some  higher  values 
of  q,  the  right  side  of  the  basic  mesh  relation  has  some  coefficients 
which  are  irrational  and  mutually  independent;  e.g.,  for  q  =  5,  after 
manipulation 


L  ■ M  =  L'2  V 3,851 852  +  853  Sm> 


■2c(,42  ■  3<S51  S52  +  s53  S54> 


(A-21 ) 


-3(S,.  S  +  S  S  ) 
51  52  53  d4 


whe  re 


c  =  cos  (2i;/10)  =  cos  36  deg 
r-  a 

y  = )  s .  s 

L  ,  L  qi  q,i+J 


(A-22) 


Since  the  lowest  order  polynomial  equation  with  rational  coefficients 
satisfied  by  cos  36  deg  is 


3  2 

8c  -  8c  +1=0 


(A-23) 


c  and  c  are  irrational  and  incommensurable,  and  each  of  the  last 
two  terms  in  the  mesh  relation  must  be  zero.  Hence 


l  -  m =y  -) 

<-50  Hi 

Hi  H2  5 


=  3tS51  S52  +  SS3  S54) 


(A-24) 


(A-25) 


Using  the  not;. ;  ion  given  by  (A-22),  the  mesh  relations  for  small  vai 
o  {'  q  a  re 


q  =  1: 


r Eq .  (32)] 


L  +  M(L  -  1) 


•I 

Ho 


q  =  2: 


q  =  3-  L  - 


q  s  4 


q  =  5: 


•M*lrL 
H<rL 
-«-y  -i 

Ho  H2 

L-L- 

t  sr  r  r~ 

q  =  6:  L  -  M  =  ^  +  )  -  >  ~  ) 

Ho  Hi  H2  H 

q  =  8:  L-M--V  -  V 

Ho  H4 

^81  ^83 


3  (S51  S52  +  S53  S54) 


q  =  12:  L 


e.g.  [Eq. 


e.g.  [Eq. 


M=y  +  y  -y  -s 

H.2,0  H.2,2  H.2,4  H.2 ,6 

l  -I 

H2 , 1  H.2  ,5 


Since  no  choice  of  signs  or  indices  has  been  made  here,  all  of 
these  re^tions  must  be  obeyed  cyclically  in  the  choice  of  se¬ 
quence  origin.  The  same  is  true  of  the  universal  auxiliary  re¬ 
lation  , 


Iv-I 

t=l  k=l 


=  ,‘-r  -  N  =  (L  +  M(L-l)) 


1/2 


(A-6)] 


(A-10)  ] 


(A-26) 
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APPENDIX  B.  THE  SHIFT-AND-ADD  RELATIONS 

If  an  m-sequence  is  shifted  t  digits  and  added  to  the  un¬ 
shifted  sequence  mod  2,  the  result  is  the  same  sequence  shifted 
to  r '  digits.  It  is  necessary  to  know  the  dependence  of  r'  on 
r  when  designing  a  linmax-coded  radar  using  the  transmitter  leak 
age  for  the  local  oscillator  (ref  25),  and  for  certain  other  pur¬ 
poses.  Here  we  state  and  illustrate  s  ime  relations  between  t  and 
r /  without  proof;  the  proofs  all  proceed  straightforwardly  using 
the  matrix  theory  found;  e.g.,  in  reference  IT. 

If  a  shift  of  T  digits  in  the  addend  sequence  leads  to  a 
shift  of  t '  in  the  sum  sequence,  we  write  r  — ►  t ' 


(B  1) 


Hence  the  relation  (  — *►  )  is  reflexive  and  should  be  written 

As  a  result  we  will  symbolize  t  and  t  more  symmetrically  by  and 


(B-2) 


■  2t2,  both  sides  modulo  L  =  2  -1. 


(B-3) 


( — t  )  -<->.(  t,,  -  r^)  ,  both  mod  L. 


If  (Tj,,  Tg)  describes  a  two-tap  linmax  SRG,  then 


for  the  m-sequence  generated  by  that  SRG. 

There  is  one  oiher  relatiot  which  permits  writing  the  (r, r  ) 
dependence  compactly: 


(B-4) 


T-  “  T,  -4-*- 7 
2  l  ;i 


<T,  '  V 


then 


(B-5) 


T4  -(T3  -(T2 


T1)B 


(mod  L) 


As  a  result  the  entire  dependence  given  in  section  4.4.2  for  the  (3,5) 
m-sequence  can  be  written  as  a  5  x  3  matrix  (generally  a 


[ 


2n  -  2 


]  x  3  matrix,  where  [  ]  symbolizes  "smallest  integer  not  less 


than.") 


where,  for  example, 


1 

2 

4 

7 

8 

1 

13 

17 


13 

26 

21 

9 

11 


17 

3 

6 

15 

12 


1  +  13  =  14 
13  +  17  =  30 
17  +  1  =  18 


and  similarly  for  each  row.  Every  integer  from  1  to  30  appears  ex¬ 
actly  once  either  as  an  element  or  the  sum  of  two  elements  in  some  row, 
the  rows  double  as  in  relation  <B-2>  and  the  sue;  of  each  row  is  31. 


The  entire  matrix  can  be  written  using  (3-4),  which  in  this  case 


reads 


followed  by  alternate  applications  of  relations  (B-2)  and  (B-3) .  This 
is  not  true  for  larger  n:  the  results  of  not  more  than  (6n)  shifts  can 
be  obtained  this  way.  Here  this  gives  ail  JO.  For  larger  n  or  for  SRG' 
with  more  than  two  taps,  it  is  necessaty  to  calculate  the  results  of 
some  shifts  by  actual  addition  of  sequences.  Relations  (B-2)  anti  (B-3) 
still  greatly  reduce  the  labor  involved,  but  by  a  factor  son  .  hat  less 
than  (6n),  due  to  short  cycles. 


APPENDIX  C:  FILTER  INTEGRATION  OF  M-SKQUENCES 


We  will  treat  a  simple  case:  the  RC  integration  of  a  mixer 
output  when  the  mixer  inputs  are  both  m-sequence  telegraph  signals. 


M-SEQUENCE 
TELEGRAPH  SIGNAL 


Figure  C-l.  Filter  integration  of  an  m-sequence. 


As  a  function  of  the  history  of  the  input  voltage  to  the 

integrator  E.  ,  tne  output  voltage  E  ^  is 
in’  out 


E  .  (t) 
<  it 


■‘f 


eK(u-t}  g  }  du 
in 


where 


-oo 

K  =  1/RC 


(' 


When  t  =  0,  E  is  dc  and  E  asymptotically  approaches  E  .  When 

t  is  a  non-zero  integer,  E  °is  the  same  m-sequence  according  to  the 

shi  ft-and-add  property.  I^nthe  sequence  is  ( - +-+-r)  and  the  recent 

historv  of  E  is  (•• - +-++ j — )  then  the  integrand  is  as  shown  in 

i  n 

figure  C-2 . 


INTEGRAND  ^"1 


At  the  end  of  any  bit,  the  integral  is 


E  =  (1  -  e"Kt°)  )  a 
out  L  n 


nKt 
e  o 


..  -Kt0.  (  -LKtQ  -2LKt* 

=  ( 1  -  e  u)  a  1  +  e  °  +  e  o  + 

o  V 


/  -Kt  - ( L+l )Kt0  ] 

+  a-l  (  6  °  +  e  °  +  ■■■  ) 


-(L-l)Kt, 


+  "-L+l  6 


(C-2) 


l-e-Kt°  fa 

l-e~LKt°  n4  &n 


nKt0 
e  ° 


(C-3) 


Unlike  the  output  of  an  L-bit  block  integrator  (one  which  cnm- 
L 

putes  )  a,  a,  exactly)  with  the  same  input,  this  E  is  not  con- 
k  k+T  out 

1 

stant  in  time.  If  the  "recent"  part  of  the  sum  is  a  section  of  the 
sequence  containing  predominantly  (+)'s,  the  sum  will  be  positive, 
etc.  Expanding  the  exponential, 

-U1  -L+l 


( 1  +  nKt  +  • • •  ) 
o 


(C-4) 


0  0 

jind  the  first  term  does  not  fluctuate;  it  is  just 

-L+l 


( C-5 ) 


A  crude  limit  on  the  fluctuating  second  term  can  be  made  by 
assuming  that  the  "recent"  half  of  the  sum  has  all  a  equal  to  (-1) 
and  the  distant  half  (+1). 


-L+l  -jL/2 

l  %  "  <  -  Z  n  +  l  " 

0  0  -L/2 


( C-6 ) 


^  La/4. 


Doubling  this  quantity  gives  a  pessimistic  estimate  of  the  peak-to- 
peak  fluctuations  of  the  sum  in  Eq  as  the  recent"  part  of  the  se¬ 
quence  changes  back  and  forth  between  "predominantly  (  +  ) 's  and  "pre¬ 
dominantly  (-)'s"as  time  goes  by. 


A  much  better  estimate  can  be  made  by  using  a  conjectured  approx¬ 
imate  limit  on  the  truncated  autocorrelation  function  of  an  m-sequence 


Then 


where 


Hence 


a.  a,  <  ( 3L/2 ) 
k  k+o  _ 


I 

t=l 


l  \n  £  £  "  -  l  " 


1/2) 


J  <  L 


1-d  1 

1/2 


d  =  (3L/2) 


y  „  a  < 


L  ( 3L/2 ) 


1/2 


( C-7 ) 


(  C-c 


( C-9  ) 


For  Perron  sequences  the  limit  on  the  truncated  autocorrelation 
seems  to  be  much  smaller. 


